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The ACS Virgo Cluster Survey V: SBF Calibration for Giant and 
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ABSTRACT 

As part of the Advanced Camera for Survey (ACS) Virgo Cluster Survey, we 
have measured Surface Brightness Fluctuations (SBF) in a sample of 100 early- 
type Virgo galaxies. Distances derived from these measurements are needed to 
explore the three-dimensional structure of the Virgo Cluster, study the intrin- 
sic parameters of globular clusters associated with the program galaxies, and 
compare with the galaxy distances derived from globular cluster luminosity func- 
tions. Our SBF measurements have been performed in the F850LP bandpass of 
the Wide Field Channel of the ACS on the Hubble Space Telescope. These are 
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the first measurements of this kind, and we present the first SBF calibration for 
this bandpass. The measured fluctuations depend on galaxy stellar population 
variations, which we quantify by galaxy color (g 475 — -2s5o)o; where g 475 and z 850 
are the galaxy magnitudes, respectively, in the F475W and F850LP ACS filters. 
We derive the following calibration for the absolute SBF magnitude M 850 : 

M 850 = -2.06 ± 0.04 + (2.0 ± 0.2) x \{g A7b - z 850 ) - 1.3] (1) 

in the range 1.3 < (#475 — ^850)0 < 1-6, and 

M 850 = -2.06 ± 0.04 + (0.9 ± 0.2) x [(5475 - ^50)0 - 1-3] (2) 

in the range 1.0 < (#475 — ^85o)o < 1-3. The quoted zero-point uncertainty here 
includes all sources of internal error; there is an additional systematic uncertainty 
of ~0.15 mag, due to the uncertainty in the distance scale calibration. 

Physically, the two different color regimes correspond to different galaxy 
types: giant ellipticals and SOs at the red end, and early-type dwarfs at the 
blue end. For the first time in SBF studies, we are able to provide a firm em- 
pirical calibration of SBF in early-type dwarf galaxies. Our results agree with 
stellar population model predictions from Bruzual & Chariot (2003) in the range 
1.3 < (#475 — z 850 ) < 1-6, while our empirical slope is somewhat steeper than the 
theoretical prediction in the range 0.9 < ((7475 — Z85o)o < 1.3. 

Subject headings: galaxies: distances and redshifts — galaxies: dwarf — galaxies: 
elliptical and lenticular, cD — galaxies: clusters: individual (Virgo) 

1. Introduction 

The Advanced Camera for Surveys (ACS; Ford et al. 1998) Virgo Cluster Survey is a 
project based on observations with the ACS on the Hubble Space Telescope (HST), aimed at 
the exploration of the properties of 100 early-type galaxies in the Virgo Cluster, the study of 
their globular cluster systems, and the reconstruction of Virgo's three dimensional structure. 
To measure the galaxy distances needed to realize these goals, we have used measurements 
of Surface Brightness Fluctuations (SBF). The SBF method was first proposed by Tonry 
& Schneider (1988) and it is based on the fact that the variance of normalized poissonian 
fluctuations of the galaxy stellar population depends on galaxy distance (for reviews see 
Jacoby et al. 1992 and Blakeslee et al. 1999). At present, SBF measurements have been 
used to determine galaxy distances out to fa 7000 km s _1 , using data from ground-based 
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telescopes and HST (Ajhar et al. 1997, 2001; Tonry et al. 1997, 2001; Jensen et al. 1999, 
2003; Blakeslee et al. 2001, 2002; Mei et al. 2001, 2003; Liu et al. 2001, 2002; Mieske et al. 
2003; and references in Mei et al. 2005, hereafter Paper IV). 

For the Virgo Cluster in particular, SBF measurements have been published by several 
groups using data from both ground-based telescopes (Tonry et al. 1990, 2000, 2001; Pahre 
& Mould 1994; Jensen et al. 1998; Jerjen et al. 2004) and HST (Ajhar et al. 1997; Neilsen & 
Tsvetanov 2000; Jensen et al. 2003). West & Blakeslee (2000) have used published ground- 
based SBF data in a first examination of the three-dimensional structure of this cluster, 
particularly the cluster's "principal axis" and its relations to the surrounding large-scale 
structure. Our new ACS dataset (for details, see Cote et al. 2004; Paper I) constitutes the 
largest, most complete, and most homogeneous SBF sample of Virgo galaxies yet assembled. 

Since the absolute magnitude of the SBF varies as a function of the stellar population 
age and metallicity, SBF distance measurements in a given observational bandpass must be 
carefully calibrated in terms of stellar population observables, typically the galaxy color. 
Using an extensive ground-based /-band SBF survey, Tonry et al. (1997, 2001) have estab- 
lished that it is possible to calibrate the SBF method using the integrated galaxy colors. 
In particular, they find that for bright ellipticals: 1) the absolute SBF magnitude Mj is a 
linear function of the galaxy (V— I)o color over the range 1 < (V—I)o < 1.3 with an internal 
scatter < 0.1 mag; 2) the zero-point of this relation is universal over this color range; 3) the 
slope of the relation is consistent with theoretical predictions from recent stellar population 
models. For observational bands other than the /-band, the relation between SBF absolute 
magnitude and galaxy colors has to be calibrated in that particular band. Ajhar et al. (1997) 
and Jensen et al. (2003) have calibrated SBF absolute magnitudes as a linear function of 
galaxy color for use, respectively, with the Wide Field and Planetary Camera-2 (WFPC-2) 
and NICMOS on the HST 

This paper provides the first calibration for SBF using the F850LP filter with the ACS, 
i.e, calibration of the absolute F850LP SBF magnitude M 850 , as a function of dereddened 
color in the F475W and F850LP filters. Observations and SBF reductions have already been 
described in Paper I, Jordan et al. (2004; Paper II) and Paper IV. In § 2 the main steps of 
the SBF reduction are briefly summarized. In § 3 the calibration method is discussed and 
the results are presented and analyzed. In § 4 the theoretical predictions from Bruzual & 
Chariot (2003) stellar population models are discussed. We summarize and conclude in § 5. 
The calibration presented here will be used in future papers to derive SBF distances for the 
galaxies in our sample and to explore the three-dimensional structure of the Virgo Cluster. 



-4- 



2. 



SBF measurements 



2.1. 



Observations and SBF reduction 



Our observations and data reduction procedures have been discussed in Paper I and 
Paper II, respectively. In Paper IV, we have addressed the specific issues related to SBF 
measurements with the ACS. Since the ACS camera presents significant geometrical dis- 
tortions, the correction of these distortions involves image resampling by an interpolation 
kernel. For our SBF measurements we use images geometrically corrected with the drizzle 
software (Fruchter & Hook 2002) using the Lancsoz3 interpolation kernel. We have shown 
in Paper IV that, when using this kernel, the distortion correction does not add significant 
error to the SBF distance measurements. The main steps of SBF measurements are also 
described in Paper IV. We summarize here the main points. 

We measure SBF in the F850LP bandpass of the ACS Wide Field Camera. The process 
involves fitting a smooth model to each galaxy in both the F850LP and the F475W band 
(Papers II and IV) and subtracting the models from the original images to produce "residual 
images" . The residual images in each band are then used to construct catalogs of external 
sources, i.e. globular clusters and background galaxies. From these catalogs we create 
a 'source mask' to remove the contamination from the sources to the surface brightness 
variance. All sources brighter than a completeness magnitude m cut , that corresponds to a 
completeness of 90% at a given radius, are masked from the image. The power spectrum of 
the residuals is then calculated in multiple concentric annuli centered on the galaxy center, 
by multiplying successive annular mask functions with the 'source mask' (i.e., the total 
mask function). The image power spectrum is the sum of two components: a flat, white 
noise power spectrum and the combined power spectrum of the SBF and undetected external 
sources, both convolved by the PSF (Point Spread Function) in the spatial domain. In the 
Fourier domain, this sum can be written as: 



where the "expectation power spectrum" E{k) is the convolution of the power spectra of 
the normalized PSF and of the total mask function of the annular region being analyzed. 
In Paper IV we have shown that this power spectrum is not significantly modified by pixel 
correlation induced by our interpolation kernel (see Paper IV for details). Armed with a 
clearer understanding of the systematic errors arising from the data interpolation, we have 
measured unbiased values of P and Pi for multiple annuli in each galaxy as described in 
Paper IV. As the minimum innermost radius, we have taken pa 1.5". However when dust is 
present in the galaxy center, the dust regions are masked (see Paper I), and the minimum 
radius without dust contamination has been taken as the minimum annular radius. For the 



P(k) = P x E(k) + Pi 



(3) 
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maximum outermost annular radius, we have chosen the radius where the galaxy brightness 
falls to half that of the sky (the range of maximum outermost radii being between m 10" 
and ~ 40"). The next step is then to extract the SBF signal from P by subtracting off the 
residual variance P r from faint sources including undetected globular clusters and distant 
background galaxies. Because of the exquisite resolution and depth of our ACS images, these 
correction are very small (see the following section). 



2.2. Contribution to the SBF from external sources 

The contribution to the power spectrum from these external sources is indistinguishable 
from the SBF contribution, since both are spatially convolved with the PSF. Following Tonry 
et al. (1990) and Blakeslee & Tonry (1995), we can write 

P = PSBF + Pr (4) 

where 

2 2 

= a gc + % ^ , g v 

(galaxy) 

o 2 gc is the contribution to the fluctuations due to globular clusters, and a\ is the contribution 
due to background galaxies. The mean galaxy brightness in the denominator is chosen for 
compatibility with the Po definition above. 

To estimate the P r contribution, sources were extracted from the same catalog used for 
masking external sources from the residual images (see above and Paper IV). Source counts 
were fitted to a composite globular cluster and external source luminosity function using a 
maximum likelihood approach. For the former, we adopted a Gaussian luminosity function, 
where N gc is a function of both radius r (distance from the galaxy center) and magnitude 
m: 

N 9C (m, r) = -j=- exp ^ — j . (6) 

For the Virgo Cluster in the F850LP filter, we expect m 9 p c eak 23 mag and a 1.35 
(Harris 1991; Whitmore et al. 1995; Ferrarese et al. 2000; Ferrarese et al. 2003? ). We use 
these values as our initial estimates but iterate on m^ ai to improve the fitted luminosity 
functions. 

For the background galaxies, we assumed a power-law luminosity function: 

N bg (m) = iV fe ° 1 7m (7) 
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Fig. 1. — The fit of total external source counts as a function of magnitude for NGC 4472 
(VCC 1226), shown as a continuous line. The globular cluster luminosity function fit is 
shown by the dashed line, and the background galaxy luminosity function by the dotted 
line. Different plots correspond to different annuli, whose range is shown in the top of 
each plot. The diamonds are the total number of external sources in magnitude intervals of 
0.4 mag. 
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with 7 = 0.35. From the ACS faint galaxy counts of Benitez et al. (2004) 7 = 0.33 in the 
F814W filter. At present, a value for 7 in F850LP is unavailable in the literature. However, 
the F814W and F850LP bandpasses are quite similar so we expect no strong differences in 
7 between these filters. A fit of faint galaxy counts in our ACS blank fields gives us a slope 
of 7 = 0.37 (Peng et al. 2005). If we vary 7 in the range 0.33 to 0.37, our fits produce 
differences of only a few thousandths of a magnitude in the final SBF magnitudes. 

In fitting the total source counts (N ogc and N b g ), we kept as fixed parameters a and 
7, and iterated on Tn 9 ^ eak , as a function of the galaxy distance. Sources identified as bright 
foreground stars were removed from the images but not included in the fits. From the 
estimated N ogc and N b g per pixel, P r was calculated as the sum of (Blakeslee & Tonry 1995) 

a 2 gc = I]V osc l0 a8[mi -<- +a4CT2ln(10)] (8) 

, m cut -m^ afc + 0.8a 2 ln(10) 

x erfc — = 

y/2a 



and 

**» ~ (0.8 - 7) ln(10) ' ( j 

where mi is the photometric zero-point on the AB system (24.862 mag for the F850LP; 
Sirianni et al. 2005) plus reddening corrections (Schlegel et al. 1998) as given in Paper II. 
An example of the external source luminosity function is shown in Fig 1 for NGC 4472 
(VCC 1226). Finally, the z-band SBF magnitudes are given by 

m 850 = -2.5 log(P - Pr) + mi , (10) 



The uncertainty on m 8 5o i n each galaxy has been calculated by propagating the measure- 
ment errors in Eq. 10 for the different annuli. This means that we have added in quadrature 
the errors due to the uncertainty in the fit of P , the uncertainty in the galaxy flux, and 
the uncertainty in P r (ss 25% P r , calculated as the variance in P r when different initial fit 
parameter are chosen). We neglect the ~ 1% systematic uncertainty in rri\ (Sirianni et al. 
2005) for the calibration, since it is common to all galaxies. The total uncertainty on the 
average m 85 o for each galaxy has been calculated as the quadrature sum of the errors for the 
different annuli which we are averaging, divided by the number of annuli. 

Finally, since m 8 5o — M 850 = 5log(d), where d is the galaxy distance, two more sources of 
uncertainty must be added in quadrature to the measurement error in m 850 for the purpose 
of this calibration. First, we include an expected "cosmic scatter" of 0.05 mag in M 850 (see 
Tonry et al. 1997), representing the intrinsic dispersion in M 850 for galaxies at the same 
observed color. Second, an uncertainty of 0.11 mag is introduced to account for the actual 
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distance dispersion from depth effects in the Virgo Cluster. This value has been derived 
from the rms scatter of ~ 3° in position on the sky of the 100 galaxies in our sample. 

3. Calibration of the ACS M 850 

3.1. Color measurements 

We seek to calibrate M 850 as a function of the galaxy color (F475W— F850LP), hereafter 
referred to (#475— z S50 ), since the ACS F475W approximates the SDSS g and F850LP the 
SDSS z. The galaxy colors are measured in the same annuli, and with the same masks, as 
those used for measuring the SBF. The sky values in the two filters for our final calibration 
sample are the same as in Paper II. Galaxies for which the sky estimates were especially 
difficult (due to strong gradients from bright neighboring galaxies) have been excluded for 
the purposes of this calibration (see below). The errors on the colors have been calculated by 
summing in quadrature the error in galaxy flux and the error due to sky subtraction. When 
the colors have been averaged on more annuli, the errors in each one have been added in 
quadrature, squared and divided by number of annuli. An additional uncertainty of 0.01 mag 
in color is added in quadrature as an estimate of the error due to imperfect flatfielding (see, 
e.g., Sirianni et al. 2005). 

The colors have been dereddened as explained in Paper II, with the reddening for each 
galaxy position from Schlegel et al. (1998). The mean reddening for the sample is E(B—V) = 
0.029 mag, with a standard deviation of 0.008 mag. The adopted extinction corrections are 
appropriate for elliptical galaxies, A g = 3.634 E(B — V) and A z = 1.485 E(B — V). Hereafter, 
we use (#475— ^850 )o to refer to the reddening-corrected galaxy color. 



3.2. Calibration Slope 

For the calibration of the J-band SBF survey, Tonry et al. (1997) divided their sample 
into galaxy groups and clusters, and assumed that, since the galaxies within an individual 
group or cluster are at nearly the same distance, SBF magnitudes will simultaneously fit the 
relation: 

m I = (m° I ) + P[(V- I) -1.15] (11) 

where (mj) is the group mean value at the chosen fiducial color (V — I)o = 1.15 mag, that 
was fitted for each group to establish the universality of the SBF calibration. To obtain the 
absolute SBF magnitude M/, Tonry et al. (1997) calibrated the zero-point in Eq. 11 with 
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Fig. 2. — Our SBF galaxy sample (85 galaxies) is shown. Elliptical galaxies are represented 
by open squares, SO by open triangles, dwarf ellipticals by filled squares, and dwarfs SO by 
filled triangles. The outliers are labeled. 
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Fig. 3. — Our final calibration sample (81 galaxies, plus VCC 21; see text) is shown, together 
with the fitted relations. Elliptical galaxies are represented by open squares, SO by open 
triangles, dwarf ellipticals by filled squares, and dwarfs SO by filled triangles. The straight 
line is the fit to the full sample, while the thick continuous lines are, respectively, fits for the 
red and blue ends. At the top left, the typical SBF magnitude error (0.12 mag, including 
the uncertainty due to the Virgo cluster depth and intrinsic SBF dispersion) is shown. 
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Cepheid distances to galaxies in their groups. Thus, they obtained the following calibration, 
valid in the color range 1.0 < (V — I) < 1.3: 

Mi = (-1.74 ± 0.07) + (4.5 ± 0.25)[(V - I) - 1.15]. (12) 

The Ajhar et al. (1997) HST/WFPC2 sample spans a range in galaxy (V— 1) color between 
1.15 and 1.3 mag. Their calibration for the WFPC2 F814W bandpass is 

M 814 = (-1.73 ± 0.07) + (6.5 ± 0.7)[(V - I) - 1.15]. (13) 

The Jensen et al. (2003) HST/NICMOS sample in the F160W (1.6 /mi) filter was calibrated 
over the range 1.05 < (V — I) < 1.24 as: 

M 160 = (-4.86 ± 0.03) + (5.1 ± 0.5)[(V - J) - 1.16]. (14) 

These three samples were heavily weighted towards bright, red, early-type galaxies. 
The Ajhar et al. and Jensen et al. samples (due the small field sizes) were mainly confined 
to the central regions of bright galaxies, which tend to be quite red. Blakeslee et al. (2001) 
concluded that this fact explained the steep slope found by Ajhar et al., since the stellar 
population models predict a steepening of the SBF-color relation for very red populations. 

Calibrations of bluer samples of dwarf early-type galaxies have been proposed by Jerjen 
et al. (1998, 2000, 2003, 2004) and Mieske et al. (2003) using ground-based observations of 
10 dwarfs coupled with the predicted color-dependence from stellar population models. 
However, until now, there has not been a sufficient sample of high-quality SBF measurements 
in dwarf ellipticals to permit a dependable empirical calibration. 

Our ACS VCS sample is the most extensive and homogeneous early-type galaxy sample 
available for SBF measurements. The galaxies span a range in color of 0.8 < (#475 — £85o)o ^ 
1.6 mag, corresponding to a range in (V— I)o between 0.94 and 1.3 mag (based on the color 
transformations from Sirianni et al. 2005). We adopt here the same technique that Tonry et 
al. (1997) have used for galaxy groups, fitting the relation: 

"%o = (rn° F850LP ) + (3[(g 475 - z 850 ) - 1.3] (15) 

for the Virgo Cluster. 

Among the 100 galaxies in our full sample, 14 have been omitted from our calibration 
sample because they are edge-on disk or barred galaxies which are especially challenging for 
SBF measurements, or because of difficult sky subtraction. These galaxies include: VCC 571, 
VCC 654, VCC 685, VCC 1025, VCC 1125, VCC 1192, VCC 1199, VCC 1327, VCC 1535, 
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Fig. 4. — Plot of our calibration sample in different color ranges (g 475 — z 850 ) cut < (c/475 — 
^85o)o < 1-6. In the top panel, the slope obtained by the fit in these different ranges is shown 
as a function of (#475 — z 850 ) cut . In the bottom panel, we show the corresponding reduced 
X 2 for the fits. Red galaxies tend to have steeper slopes 2) than the slope fitted to the 
overall sample (1.3 ± 0.1). 
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VCC 1627, VCC 1720, VCC 1857, VCC 2048, and VCC 2095. One galaxy in the sample, 
VCC 1499, has (#475— z 850 ) < 0.8, and is too blue to be reliably included in a calibration. 
There are also several outliers that do not appear to lie on the average SBF magnitude 
versus color relation for the cluster. Fig. 2 shows the SBF and color measurements for the 
remaining 85 galaxies. 

To calibrate our SBF measurements in the F850LP filter as a function of the average 
galaxy (#475 — z 850 ) color, we calculate average SBF measurements for each galaxy and fit 
the relation in Eq. 15 to average galaxy colors, as Tonry et al. (1997) have done for each 
group in their sample. All galaxies included in the fit have an uncertainty on fnsm less than 
0.15 mag (not considering the uncertainty due to the 'cosmic' scatter on the SBF calibration; 
see above). Errors on average galaxy colors are typically ~0.01 mag. Both uncertainties on 
magnitude and color are included in the x 2 calculations below. 

To exclude the outliers from the final calibration, we first make a linear fit to Eq. 15 
to the overall sample with a robust least absolute deviation technique (Press et al. 1992). 
We obtain (m° FS50LP ) = 29.09 ± 0.02 mag and /? = 1.25 ± 0.10. We identify four outlier 
galaxies from iterative 3-a clipping. The four 3-a outliers are: VCC 21, VCC 538, VCC 575, 
and VCC 731. Of these, VCC 21 contains numerous bright young star clusters, indicative 
of recent star formation, while VCC 731 (NGC 4365) is known from the ground-based SBF 
survey to be a Virgo member which lies on the far side of the cluster; we suspect that 
the other two may be similar cases. After excluding these four outliers and the very blue 
VCC 1499, our final sample consists of 82 galaxies. 

In Fig. 3 our final galaxy sample and an overall fit (the continuous straight line) are 
shown. Both uncertainties on magnitude and color have been included in the fit x 2 estima- 
tion. We obtain: 

m 850 = (29.09 ± 0.02) + (1.3 ± 0.1) x [(^475 - ^50)0 - 1-3] (16) 

in the overall color range 0.9 < (c/475 — ^850)0 < 1-6. This color range includes both bright and 
faint galaxies. Various lines of evidence have suggested that the slope of the /-band SBF- 
color relation is steeper for redder galaxies: for instance, the steep empirical slope found 
by Ajhar et al. (1997) for the red centers of bright ellipticals, which was also reproduced 
by the reddest composite stellar population models of Blakeslee et al. (2001). Conversely, 
the few dwarf galaxies in the Tonry et al. (1997) sample indicated a significantly shallower 
slope at the blue end, and the SBF measurements for Galactic globular clusters (Ajhar & 
Tonry 1994) showed no dependence on color. Most recent models also show a flattening of 
the relation for bluer populations. 

Our sample, more clearly than any previous study, suggests that changes in color range 
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Fig. 5. — Best-fit slopes (top panel) and reduced x 2 values (lower panel) in two different color 
ranges [blue: 0.9 < (#475 - ^850)0 < (#475 - z 850 ) div and red: (g 475 - z 8 5o)div < (#475 - ^850)0 < 
1.6] are shown as a function of ((7475 — zs5o)div The dot-dashed and dashed lines correspond 
to the red and blue sample, respectively. The solid line in the lower panel shows the overall 
reduced x 2 f° r the full sample. Red galaxies (mostly giants) prefer a slope ~ 2, while the 
blue galaxies (mostly dwarfs) prefer a slope ps 1. 
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affect the fitted slope, so we have studied the stability of the above relation as a function of the 
lower color limit, ((7475 — z 850 ) cut . When we repeat the fit in different ranges (#475 — z 850 ) cut < 
(#475 — -285o)o < 1-6, we obtain the results shown in Fig 4. On the top of the figure the fitted 
slope (3 is shown as a function of (#475 — z 850 ) cut . On the bottom is shown the corresponding 
value of x 2 - Red galaxies tend to have slopes that are steeper 2) than those obtained by 
fitting to the overall sample. This suggests that our sample would be better represented by 
two distinct samples of red and blue galaxies, each one having its own calibration slope. 

We fit Eq. 15 in two complementary color ranges : 0.9 < ((7475 — ^85o)o < (#475 _ Z85o)div 
(blue) and (#475 - z 850 ) div < (g 475 - z 850 ) < 1.6 (red), with 1.15 < (#475 - z 850 ) div < 1.35 (as 
suggested from Fig. 4). In Fig. 5 the fitted slopes (top) and the respective x 2 (bottom) are 
shown as a function of ((7475 — z 850 ) div . Dot-dashed lines correspond to the fit in the color 
range (#475 — z 8 5o)div < (#475 — ^850)0 < 1-6 and dashed lines to the fit in the color range 
0.9 < (#475 — ^50)0 < (#475 — z 8 5o)div The continuous line is the total x 2 from both fits. The 
minimum total x 2 is obtained at ((7475 — z 850 )di V = 1.26, where the red slope is 1.95 ± 0.24 
and the blue slope is 0.90 ± 0.23. The two slopes cross at (#475 — z 850 )di V ~ 1-3, which we 
take to be the color at which the two regimes separate. The x 2 f° r the total fit is minimized 
for both slope and zero-point, and the errors in both magnitude and color have been taken 
into account. 

As in Tonry et al. (1997), we find the SBF in SO and ellipticals to follow the same 
trend as a function of (#475 — ^85o)o- Interestingly, as Figs. 2 and 5 show, the dividing color 
°f (#475 — <285o)o — 1-3 is close to that which separates the two most basic galaxy types in 
our survey: giant and dwarf early-type galaxies. In fact, in our final sample, of 37 galaxies 
with color (g 475 — z 850 ) > 1-3, only one is classified as dwarf. For 1.2 < (#475 — -2350)0 < 1-3, 
10 galaxies are dwarfs and 8 are giants. For (g 475 — z 850 ) < 1.2, of 26 galaxies, 20 are 
dwarfs. Physically then, the two different calibration slopes correspond to different galaxy 
types (giants corresponding to the red end, and dwarfs to the blue end). 

Jerjen et al (2000, 2004) adopted a theoretical calibration to dwarf elliptical galaxies, 
based on a two branch prediction in the R-band SBF versus the (B — R) galaxy color from 
Worthey et al. (1994) theoretical stellar population models using the Padova isochrones. 
While the scatter becomes larger for the bluest dwarfs, we see no evidence for distinct 
branches in our large data set; nor do the Bruzual & Chariot (2003) models predict such 
branches in z 850 vs (#475 - z 850 ) . 

Our sample is sufficiently complete at the blue end that we are able to provide, for the 
first time in SBF studies, an empirical calibration of SBF in early-type dwarf galaxies. 

Splitting the galaxy sample at ((7475— Z85o)o = 1-3, we derive the following calibration 
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Fig. 6. — Average and standard deviation of m 8 5o in 0.05 mag color bins (top) and the color 
histogram (bottom) of our final calibration sample (81 galaxies, plus VCC 21; see text) are 
shown. 
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for m 850 : 

m 850 = 29.03 ± 0.03 + (2.0 ± 0.2) x [(g 475 - z 850 ) - 1.3] (17) 
in the range 1.3 < (g i75 — z 850 ) < 1.6, and 

m 850 = 29.03 ± 0.03 + (0.9 ± 0.2) x [(#475 - ^850)0 - 1.3] (18) 

in the range 0.9 < (#475 — ^50)0 < 1-3. The reduced x 2 values for the red and blue subsam- 
ples are 0.99 and 0.8, respectively. The value of the total reduced x 2 is calculated as the 
appropriately weighted sum of the x 2 values from the two independent fits. When normal- 
ized by the total number of objects used in the fits, its value is 0.9. The reduced x 2 for the 
blue subsample is lower than that for the red sample. This might be due to the difficulty in 
accounting for the intrinsic distance dispersion within Virgo, or to overly conservative error 
estimates for the dwarfs. One of the outliers rejected above, VCC 21, is no longer an outlier 
with respect to the final, shallower, blue-end calibration, and we have added it to the final 
calibration above, by iteration. 

We show in Fig. 6 the average and standard deviation of mgso in ((7475 — £ 8 5o)o color 
bins (top) and the color histogram (bottom) of our final calibration sample (82 galaxies, plus 
VCC 21). From the scatter in m 8 5o, we are aware that our calibration has a large scatter for 
colors (#475 — z 8 5o)o < 1- We will limit our absolute magnitude calibration to the color range 
1.0 < (#475 - z 85 o)o < 1-6. 

We observe a similar behavior of m 8 5o as a function of the (#475 — z 850 ) color for the 
multiple annuli within individual galaxies. However, the precise value of the slope is less 
robust and more sensitive to errors in the sky. 

3.3. Calibration Zero-Point 

Empirical calibration of the zero-point of the SBF relation requires computing the SBF 
absolute magnitude M^ 850 lp from independent distance measurements. The ideal situation 
would be to have Cepheid distances for multiple galaxies in our sample, however there are no 
galaxies in our sample with Cepheid distances. We thus calibrate our zero-points to the Virgo 
Cluster distance modulus found in the /-band SBF survey of Tonry et al. (2001), but revised 
in accordance with the final set of HST Key Project Cepheid distances from Freedman et al. 
(2001). Tonry et al. (2001; Table 4) reported 31.15 ± 0.03 mag as the Virgo SBF distance 
modulus, where the SBF zero-point was calibrated from SBF measurements for individual 
early-type spiral galaxies with Cepheid distances given by Ferrarese et al. (2000). 

However, using the final set of Key Project distances from Freedman et al. (2001), 
who adopted the Udalski et al. (1995) period-luminosity relation and added a metallicity 
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correction, the Tonry et al. (2001) distance modulus for Virgo has to be corrected with a 
shift of —0.06 mag, and becomes 31.09 ± 0.03 mag (see discussions by Ajhar et al. 2001; 
Blakeslee et al. 2002). Here, the errorbar reflects the internal error on the mean Virgo SBF 
distance; the uncertainties in tying the ground-based SBF distances to the Cepheid distance 
scale, and in the zero-point of the Cepheid distance scale itself, contribute another ~0.15 
mag systematic error. If the metallicity corrections to the Cepheid distances tabulated by 
Freedman et al. are not applied, then the Tonry et al. (2001) Virgo distance modulus would 
shift instead by -0.16 mag to 30.99 ± 0.03 mag (see Jensen et al. 2003). For our SBF 
calibration, we will adopt 31.09 ± 0.03 mag as the Virgo distance modulus. 

Thus, scaling the zero-points from Eq.17, we obtain 

M 850 = -2.06 ± 0.04 + (2.0 ± 0.2) x [(g 475 - z S50 ) - 1.3] (19) 
in the range 1.3 < ((7475 — z$5o)o < 1-6, and 

M 850 = -2.06 ± 0.04 + (0.9 ± 0.2) x [(g 475 - z S50 ) - 1.3] (20) 
in the range 1.0 < (# 475 - z 850 ) < 1.3. 



4. Predictions from stellar population models 

We now compare our results to predictions from the Bruzual & Chariot (2003) stellar 
population models, which we have calculated from the models and stellar luminosity func- 
tions kindly provided by S. Chariot. In Fig. 7 we show our data as compared to Bruzual & 
Chariot (2003) model predictions for different stellar populations. The lines connect models 
of the same metallicity: the red line is for solar metallicity, the green for 40% solar, the blue 
for 20% solar, the light blue for 2% solar, and the orange for 2.5 times solar metallicity. The 
dots on the lines represent, respectively from left to right, ages of 1, 3, 5, 8, 12 and 15 Gyr. 

The Bruzual & Chariot models predict M 850 . To compare these predictions with the 
data, we add to the theoretical M 850 a distance modulus such that the predicted average 
m 850 for solar and 40% solar metallicity at (g 475 — ^50)0 = 1-2 mag is equal to the average 
observed m 850 at the same (#475 — z 850 ) f° r our calibration sample. This procedure yields 
a theoretically predicted SBF Virgo Cluster distance modulus of M 850 — m 850 = 30.94 mag, 
which is entirely independent of the Cepheid distance scale, but we caution that it depends 
on the color range we have chosen for shifting the models into agreement with the data. 

This theoretically calibrated SBF Virgo Cluster distance modulus differs by ^0.15 mag 
with respect to the revised Tonry et al. (2001) Virgo distance modulus of 31.09 mag adopted 
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Fig. 7. — Theoretical predictions from Bruzual & Chariot (2003) stellar population models 
are shown. The red line is for solar metallicity, the green for 40% solar, the blue for 20% solar, 
the light blue for 2% solar, and the orange for metallicity 250% solar. The dots on the lines 
represent, respectively from left to right, ages of 1, 3, 5, 8, 12 and 15 Gyr. 
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above. However, the two are consistent within the uncertainties from the zero-point of 
the Cepheid distance scale, the anchoring of the ground-based SBF measurements to the 
Cepheids, and the various uncertainties in the stellar population models (see Blakeslee et 
al. 2001 for a discussion). It is also worth noting that if we had adopted the Virgo distance 
modulus of 30.99 ± 0.03 obtained by calibrating SBF to the Freedman et al. (2001) Cepheid 
distances that use the new period-luminosity calibration but omit the metallicity correction, 
then the difference with respect to the model-predicted distance modulus would be only 
0.05 mag. 

Our sample of early-type galaxies spans a wide range in metallicity and, presumably, 
age. Like the observations, the theoretical predictions also show a steeper slope at the 
red end, and a shallower at the blue end. The slopes predicted from linear fits to the 
overall model predictions are 0.0 ± 0.5 for 0.9 < (g 475 — z 850 ) < 1.25 mag, 2.6 ± 0.9 for 
1.25 < (#475 - z 85 o)o < 1-6 mag, and 1.26 ± 0.3 for 0.9 < (#475 - -^50)0 < 1-6 mag (in very 
good agreement with our empirical slope over the full range). If we fit only metallicities 
less than or equal to solar, at the red end the slope will be 2 ± 1, compared to 0.4 ± 0.4 
in the blue. The agreement between theoretical predictions and our empirical calibration is 
somewhat better at the red end. This is likely due to better accuracy in the modeling of 
highly evolved populations. The SBF predictions for young to intermediate-age populations 
is less certain due mainly to the uncertainties in the modeling of the asymptotic giant branch 
contributions. Furthermore, models of old stellar populations are better constrained by 
observations of Galactic globular clusters. For these reasons, we prefer our empirical blue- 
end slope, although additional SBF measurements of large homogeneous samples of dwarf 
galaxies (e.g., in the Fornax cluster) would be helpful in refining the SBF calibration for blue 
early-type galaxies. 

Considering the uncertainty in the Virgo distance modulus, we conclude that our z 850 
SBF data for the giant early-type galaxies are well reproduced by the models for stellar pop- 
ulations ages greater than 3 Gyr and metallicities in the range between 40% and 250% solar. 
The data for the dwarf galaxies are consistent with either much younger (less than 3 Gyr) 
solar metallicity populations or intermediate-to-old stellar populations with metallicities in 
the range 20%-40% solar. 

5. Conclusions 

As part of the ACS Virgo Cluster Survey, we have observed 100 early-type galaxies in 
the Virgo Cluster. We have measured SBF magnitudes in these galaxies with the aim of 
determining their distances and the three-dimensional structure of Virgo. We have presented 



- 21 - 



the first calibration of the SBF method using the ACS instrument, giving the F850W SBF 
magnitude as a function of galaxy (c/475 — z$ 50 )o color. Obvious outliers, most likely due 
to a real dispersion in distance for galaxies in the Virgo Cluster region, have been omitted 
from this calibration, leaving 82 galaxies with ((7475 — z S50 ) colors ranging between 0.9 and 
1.6 mag. Over this color range, the blue (mainly dwarf) and red (mainly giant) galaxies 
follow different linear relations, with the slope being shallower at the blue end. For the first 
time in SBF studies, our sample allows us to empirically calibrate SBF for dwarf early-type 
galaxies. 

By comparing with stellar population models from Bruzual & Chariot (2003), we find 
that the slope of the model prediction for msso with (#475 —£850)0 color is highly consistent 
with the observed relation for the red galaxy subsample. There is also general consistency 
at the blue end, though the stellar population models with young ages and low-metallicities 
predict a somewhat shallower slope than observed. This underscores the importance of em- 
pirical calibrations, as well as more realistic composite populations in any effort to model the 
behavior of the SBF method theoretically, especially for dwarf galaxies, for which the scatter 
in stellar population properties is larger. A forthcoming paper will present the distances to 
the individual galaxies derived with this calibration, and examine the issue of Virgo's three 
dimensional structure and its implications for galaxy and globular cluster properties. 
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